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We study the production of prompt photons in strongly coupled out-of-equilibrium Super Yang- 
Mills plasma using the AdS/CFT correspondence. Our goal is to determine the photon emission 
spectrum at different stages of a thermalization process, which is modeled via the gravitational 
collapse of a thin spherical shell in AdSs space. Particular emphasis is placed on the limit of large 
frequencies, which we are able to treat analytically. 



Introduction and setup. In heavy ion experiments 
highly energetic electromagnetic probes, such as pho- 
tons and dileptons, provide important information on the 
evolution and eventual thermalization of the produced 
plasma (for reviews, see e.g. [HE])- This is due in partic- 
ular to their weak coupling to the plasma constituents, 
which implies that once produced, they are free to propa- 
gate through the plasma almost unaltered. The spectrum 
of direct photons measured in central Au-Au collisions 
indeed shows an enhancement above the scaled p-p spec- 
trum for transverse momenta less than 2 GeV [3HS] . This 
appears to be consistent with the production of thermal 
photons in a strongly coupled quark-gluon plasma, as in- 
dicated by both RHIC [5] and LHC experiments [Tj. 

On the theory side, the dynamics of strongly inter- 
acting field theories out of thermal equilibrium is a no- 
toriously difhcult problem, as perturbative results are 
only applicable in the limit of asymptotically high en- 
ergies, while lattice methods are in general constrained 
to equilibrium quantities [SI [9]. In this context, the 
gauge/gravity duality has proved itself highly useful, as 
it allows one to reduce the strongly coupled dynamics 
of certain field theories to classical gravity problems in 
curved spacetime [T0Hl2] . Thermalization on the field 
theory side is argued to be dual to the formation of a 
black hole in an (asymptotically) Anti de Sitter (AdS) 
spacetime, which one can model e.g. via the gravitational 
collapse of a thin shell of matter in this geometry [I3HI8] . 

In the present paper we concentrate on the production 
of prompt photons in a thermalizing plasma, aiming to 
extend our previous work on dileptons [17]. In particular, 
we are interested in the (isotropic) differential production 
rate of on-shell photons with energy = k [THl [2D] j 

fco^ = -L^ ^ " ^M-n< (fco = k) , (1) 

d^k Airk dko ' p^^v y : v / 

where a is the fine structure constant and II^j, the elec- 
tromagnetic current Wightman function. In thermal 
equilibrium, the fiuctuation dissipation theorem allows 
one to further relate 11^^ to the (transverse) photon spec- 
tral function Xfii^^ 

ry^^n<,(fcO = fc) = nB(fc°)x;^:(fco) ^ nB{k')x{ko) ■ (2) 



While this is not generically true out of equilibrium, it 
was shown to hold in the quasistatic limit of the falling 
shell system in appendix B of [T7] . 

The gravity setup we work in is thoroughly explained 
in [17j, and only briefly summarized here. We place an 
infinitesimally thin shell of unspecified matter at some 
radius in AdSs space, and let it fall gravitationally 
in the radial direction. The limit where Vs approaches 
the Schwarzschild radius rh of the shell is conjectured to 
correspond to thermalization in strongly coupled, large- 
Nc Af = 4: Super Yang-MiUs (SYM) plasma, while rn 
(and thus the mass of the shell) is related to the final 
equilibrium temperature of the field theory. 

Next, we introduce in the SYM theory a U(l) gauge 
field coupled to the conserved current corresponding to a 
U(l) subgroup of the SU(4) R-symmetry. According to 
the standard AdS/CFT prescription [5T], the correspond- 
ing retarded correlator — and thus the 'photon' spectral 
function — is calculable by studying linearized fiuctua- 
tions of an (unrelated) U(l) gauge field in the bulk, where 
particular care must be taken of the boundary conditions 
of the field at the shell [T7]. We perform the calculation 
in the so-called quasistatic approximation ,13J, in which 
the shell is assumed to move slowly compared to the other 
time scales of interest. As discussed in section 4 of [T7] . 
this condition is satisfied in the later stages of thermal- 
ization as well as for highly energetic photons, with fre- 
quency larger than the inverse timescale associated with 
the motion of the shell. 

Finally, during the past few months, a number of re- 
lated works on the subject of holographic thermalization 
have appeared on the arxiv; for details, we refer the in- 
terested reader to refs. [^21 - BS] and references therein. 

The metric and the equation of motion. The metric 
of the falling shell setup consists of a black hole solution 
outside the shell, r > r^, and of pure AdS space inside it, 
r < Tg. The parameter is assumed to be larger than 
the Schwarzschild radius of the shell, > r^, so that 
initially there is no black hole and the dual system is out 
of equilibrium. Assuming to be much larger than the 
curvature radius of AdS space, ^ L, the line element 
of the spacetime can be written in the form 

ds'-'^{-fir)dt^ + d.^) + ^^^, (3) 



2 



where we have defined 



fir) = 



1 - 
1, 



for r > 
for r < r. 



(4) 



The Hawking temperature of the eventual black hole is 
given by r = this is also the temperature in the 

field theory that we will be referring to throughout the 
evolution of the system (although it is in principle well- 
defined only at the end). For notational convenience, we 
will in the following set L = 1 and introduce a new radial 

variable u = In this new coordinate, the boundary 
is located at u = and the horizon at u ~ 1, while the 
only relevant parameter of the gravity setup, the shell 

location, takes the form = 

To study photon production, we need to investigate 
a transverse electric field E-^{t,u,'K.) [2Ql[27]. Denoting 
ko — Lo and using the fact that for an on-shell photon k — 
fco, we obtain for its equation of motion in momentum 
space dU 



duf 



duE^ 



-E^=0 



d 



(5) 



where we have introduced the dimensionless variable uj = 
^ = This equation has two linearly independent 

sotutions, which are given by hypergeometric functions 



] (see also [20]) 

Et{u,u) = (l-u)- 
( 1 + i 



2^1 



'{l+U 

l + i 



oj;l — iuj; 



2^1 



F, 1- 



1 



i 1 

-cj,-- 



-w; 1 + iuj; 



(l-u) 



(1 



(6) 



(7) 



^out 



an outgoing 
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where E^^ satisfies an infalling and iJ^ 
boundary condition at the horizon. 

In a black hole background, the physical solution is the 
one obeying infalling boundary conditions, as classically 
nothing can escape from a black hole [51] [57] . In the case 
of a falling shell, there is, however, no horizon, and hence 
we must write our full solution as a linear combination 
of the infalling and outgoing components, 

^^utsidol'^, u) = c+E^^{lj, u) + c^E^^^iuj, u) , (8) 

and match this to the solution inside the shell. As dis- 
cussed in [T7], the discontinuity of the time coordinate 
across the shell implies that the frequency outside the 
shell is related to the inside one through a rescaling by 



V fm — 



Taking this into account, the equation 



of motion inside the shell is given by 



i7iside 



fn 



inside 



= 



(9) 



the solution of which reads in terms of Besscl functions 

[Ji (2^7^) + iYi {2iL^/u)] 

u^O 1 .r./^^ ^^^^ 




FIG. 1: The trace of the spectral function x('^j normalized 
by ljN^T'^. In order of increasing amplitude of the fluctua- 
tions, the curves correspond to Va/rh ~ 1.001 (green), 1.01 
(red) and 1.1 (blue). The dashed black line stands for the 
thermal spectral density. 



where we have denoted 




(11) 



The matching of the two solutions then leads to the result 
(for details, see [17] ) 



-^in^M ^inside V7m-E'insidc^"^in 



^+ -^rat^u ^inside VTm-Bjl^sidc^M-^iit 

In the limit ^ 1, this ratio vanishes as 

i(l-U,)l/2-« 







(12) 



(13) 



which serves as a consistency check of our calculation 
(note that we are only interested in real values of w here) . 

Photon production. From eqs. ([T])-([2]), we see that in 
order to study the production of prompt photons, it suf- 
fices to evaluate the trace of the transverse photon spec- 
tral function. As the latter is given by the imaginary part 
of the retarded R symmetry correlator n-*^ , we can follow 
the standard AdS/CFT prescription of [21] and write 



J ^" -^OTtside 



(14) 



where we have used the shorthand 0(w = 2'itTw) = 

for the functions x s^nd 11-'- . This quantity is furthermore 

conveniently expressed in terms of the Wronskian of 



-^^utside(^.") 
^^o'-ut.ide('^." = 0) 
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Imn-L(tj) = ^ W{uj,u), 



(15) 



(16) 



W{uj,u) 



Im f{u){F^iid,u)yduF^{uj,u) 



FIG. 2: The relative deviation of the photon spectral function 
from its thermal limit, R{lj,Us), for Vs/rh ~ 1.001, 1.01, 1.1. 
The color coding is as in fig. [T] 



FIG. 3; The photon emission spectrum dT^/ko, normalized 
by aN^T^, for rs/rh = 1.001, 1.01 and 1.1. The color coding 
is as in fig. [l] 



The advantage of this formulation is that it can be easily 
shown that W{uj, u) does not depend on the variable u, 
i.e. W{6j,u) = W{6j). Thus, we can determine its value 
in the numerically more convenient limit m — ?> 1 instead 
of the boundary, u = 0. 

A lengthy but straightforward calculation employing 
the formalism explained above leads us to the result 



1 



D{Qj) 



(17) 



where we have reinstated the explicit Us dependence of 
the spectral function (originating from c±) as well as de- 
noted 



D{u:) ^ \F\^[l 



2 Re 



2 



(18) 



1 + i . 1 - i . 
F = 1 1 + ^^^-^ 1 - itj; -1 



From eq. (13 1, we know that the ratio c_/c+ approaches 
when — > 1. In this limit, the spectral function thus 
reduces to its known equilibrium limit. 



Xthormal ij^) 



N^T^jj 1 



(19) 



as required by consistency (cf. eq. (3.18) of [TU]). 

The novel aspect of our result is that eq. (T7|) gives 
the photon spectral density both in and out of thermal 
equilibrium. As thermalization is parameterized by the 
shell location, we use several different values of when 
plotting this function in hg. [ij In fig. [2] we also show the 
relative deviation of the spectral density from its thermal 
limit, 



i?(w. Us) 



X{>^,Us) - Xthcrmal(t^) 
Xthermal(w) 



(20) 



which exhibits an oscillatory pattern much like the one 
already observed in the dilepton spectral function in |17| 
(for a discussion of the physical origin of the oscillations. 



see the same reference). As the shell approaches the hori- 
zon, Ts — >■ r/i, the amplitude of these oscillations is seen 
to decrease, as it of course should. The amplitude is also 
seen to be dampened at large values of di, which we in- 
terpret as a sign of the usual top/down type pattern of 
holographic thermalization. In the opposite limit of small 
Lo, the out-of-equilibrium x/'^ is on the other hand seen 
to vanish; this is, however, merely a sign of departing the 
range of validity of the quasistatic approximation. 

In fig. [3] we finally show the photon emission rate per 
unit volume as a function of w, 



dkn 



ak 



'nBiko)x{ko) 



ko=k=27TTui 



(21) 



We observe a slight enhancement of the photon produc- 
tion rate, when the system is out of thermal equilibrium. 
Nearly all traces of the oscillatory behavior have, how- 
ever, been suppressed by the Bose-Einstein distribution 
function, and in particular the overall shape of the spec- 
trum stays largely unaltered when rg > r/j. It is in- 
teresting to contrast this behavior with that of the lead- 
ing corrections to the equilibrium photon emission 
rate studied in |29| . While the equilibrium emission spec- 
trum of [19] was seen to get a moderate enhancement also 
in this case, it was found that the peak of the function 
moved towards the IR with decreasing A, while in our 
case it shifts towards the UV. 

Large Co limit. In analogy with the dilepton case of 
ref. [T7], it would clearly be of some value to get an 
analytic handle on the behavior of the photon spectral 
density in the limit of large w, where the quasistatic ap- 
proximation is known to work. To this end, we follow the 
discussion of appendix A of 1W\ and first note that for 
large w, the general solution to the equation of motion 
([5]) can be expressed in terms of Airy functions. 



C2(a;)Bi(^2/3C(-^i)) 



(22) 
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with Ci{ili) denoting unknown coefficient functions and 

2/3 



Ciy) ^ 

g{y) = i\n 

Kv) 



r/3[ 3.9(~y) 

I 4 



-y-i 



V 



(23) 

— 2i arctan(-v^— y) — tt , (24) 
(25) 



(i-m(y) ■ 

One sliould note liere that the quantity inside the curly 



brackets in eq. ( 23 ) is real and positive in the range < 
2/ < 1- 

It is a straightforward exercise to show that close to 
the horizon, u w 1, eq. (22 1 reduces to 



E. 



(^, u) ^ ^^^^ Ci(^)2^'^e— ^/4(1 ~ «)-^"/2 



-C2(w)(2i-*'^e''^"/'*(l-w)' 



-i2*"e-"^"/4(l-w)-*'^/2> 



2^/2 r _ 

___ |2-e— /4(Ci(^) - ^C2mEU^, u) 
+2'-'''e"''/^C2{io)Eout{uJ,u)} , (26) 



where in the latter equality we have used the w — > 1 limit 
of eqs. ([6|-([7]). Denoting finally 



a 



2(i/2g-i-7r/4 

2v^wi/6 



/3 = 2*'^e-""/^ , (27) 



we then obtain upon comparison with eq. ([8]) the coeffi- 
cients Ci{(jj), 



(28) 
(29) 



A very useful byproduct of the above calculation is 
that we can now write the properly normalized infalling 
and outgoing solutions to the equation of motion ([s]) in 
the large-cl; limit as 

Etiu^^u) = I^^^pgLMiu^^i-u)), (30) 



Eout{<^,u) 



+Bi(^2/3^(_^)) 



(iMico'/^Ci-u)) (31) 



from which we further obtain 

Ei, (lJ, u) "^"^ 2-'^(l±2^)/2g±»^(.+23(«))/4^-l/4 (32) 
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FIG. 4: The behavior of the full function R{(2j,Us) (blue) to- 
gether with its IsLTge-Lj limit i?asym('i, lis) (red) for r^/rh = 
1.1. 



This is helpful in determining the large-w limit of the 
retarded Green's function of interest, 



nasym('^,Ws) = - 



N^T^ 3l/3p(2/3) 



1 



; r(l/3) 

22i(ig-i7r/2([i + l/3) < 



(33) 



1 + 2'^i'^e—i'''/^{'^—'^/3) 



2/3 



which clearly reduces to the correct equilibrium result of 
[H] when c_/c+ 0. 

Inserting finally the straightforwardly obtainable large- 
w limit of c_/c_|_ to eq. (33 1, we obtain an analytic, 
though somewhat cumbersome-looking expression for the 
retarded correlator, 

N^T^ r(2/3) 



ni 



8 r(i/3) 



X 6^2-/331/3^2/3 (34) 



o 3/2 . 
3/2 



^1 - m2(1 + ^1 -u2)(i + e»-/3e*<:'s(".)) 

Sul^^cb - ^1-^2(1 + ^1 -u2)(i - e»2-/3e»'iff("»)) 

From here, we easily obtain the large-w limit of corre- 
sponding spectral density, 

^2y2 35/6p(2/3) 



r(i/3) 



yr^(i + yr^) 



8m: 



3/2 



cos{ujg{us)) + \/3sm{ojg{us)) 



(35) 



which is seen to reduce to the correct equilibrium limit as 
Us — > 1. In fig. |4] where we plot the asymptotic form of 
R, derived using eq. ( 34 1 , together with our full numerical 



result, we observe excellent agreement of the two curves 
already at moderately small frequencies. 

Conclusions. In the paper at hand, we have investi- 
gated the production of prompt photons in strongly cou- 
pled, out-of-equilibrium J\f — A SYM plasma. This work 



5 



can be viewed as a generalization of the equilibrium com- 
putation of Caron-Huot et al. [19] to a specific model of 
holographic thermalization, involving the gravitational 
collapse of a thin spherical shell in AdSs spacetime [13] ■ 

Our main results are depicted in figs. [l}{3j They dis- 
play a distinctive pattern of fluctuations in the spectral 
density, which however are significantly dampened in the 
photo-emission spectrum. In these figures, the thermal- 
ization process is parameterized by the radial location of 
the shell, but as discussed in section 4 of [17], explicit 
time dependence can be introduced at any point by solv- 
ing the equation of motion of the shell (upon specifying 
its matter content and initial condition). 

The fashion, in which the photon spectral density is 



seen to approach its thermal limit, is very similar to that 
computed for dileptons at rest in ref . [17] . It is indicative 
of a top/down type thermalization pattern, and shows no 
apparent dependence on the virtuality of the produced 
photons (see cf. e.g. ^30] for a different conclusion in the 
case of jet quenching). 
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